The Big Trip and Wheeler-DeWitt equation 
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Abstract 

Of all the possible ways to describe the behavior of 
the universe that has undergone a big trip the Wheeler- 
DeWitt equation should be the most accurate - pro- 
vided, of course, that we employ the correct formula- 
tion. In this article we start by discussing the standard 
formulation introduced by Gonzalez-Diaz and Jimenez- 
Madrid, and show that it allows for a simple yet efficient 
method of the solution's generation, which is based on 
the Moutard transformation. Next, by shedding the un- 
necessary restrictions, imposed on aforementioned stan- 
dard formulation we introduce a more general form of 
the Wheeler-DeWitt equation. One immediate predic- 
tion of this new formula is that for the universe the 
probability to emerge right after the big trip in a state 
with w — Wq will be maximal if and only if wq = — 1/3. 



Keywords phantom energy, 
tion, big trip 



1 Introduction 



Wheeler-DeWitt equa- 



In recent years the mathematical community of cos- 
mologists has provided us with a multitude of new 
amazingly original models, not the least of which 
is the "Phantom fields" concept, that (among other 
things) notably postulates the violation of the weak 
energy condition (WEC) p > 



+ p/c 2 > 
(jCaldwell et al.l (|2003h . ICaldwelll (|2002l) ). where p is the 
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flui d density and p is the p ressure (for recent reviews 



see [Bamba et al 
(12003 
(200 



( 20121) Caldwell &: Kamionkowski 



Frieman fc Turnen (|2008l ). ISilvestri fc Trodden 
Li et al.l (1201 ll)). S uch fields, as follows from 



quantum theory ICarrol et al.1 (|2003h . should manifest 



themselves via the scalar field with the negative kinetic 
term. The thorough investigation asserts that fields 
with such a property should be anything but funda- 
mental. Still, one cannot guarantee that no Lagrangian 
with the negative kinetic terms would arise to serve as 
an effective model. And that, indeed, is exacly what 
happen s in a t least some of the supergravity models 
( Nilles (|l984j ) ) . as well as in the gr avity theories with 
highest derivatives ( Pollock ( 1988)) an d even in the 
field string theory (jAref 'eva et al. in the mod- 

els describing the fermion NSR-string with regard for 
GSO_-sector (see also|Sen| (120051) 11. Lately, the "phan- 
tom energy" had even managed to penetrate into the 
domain of yet another "exotic" cos mological model - 
the brane theory (see, for example, ISahni fc Shtanov 



( 20021) . lYurov et al 



(120061) ). 

Remark 1. Although most theoreticians believe that 
the phantom energy is nothing but an effective model, 
this shared notion as of yet remains to be backed up 
by anything rigorous. In fact, there are some evi- 
dence of quite the opposite. One particular example 
is the phenomena of "crossing of the ph antom d ivisio n 



line", first dis c overed and discussed in lYuroyl (2003) 



Elizalde at al.l (|2004l) . lYurov fc Ver cshchagin 



2004) 



and later interpreted in lAndrianov et al.l ([20051 ) . This 
"crossing" is smooth and it is an inherent property 
of some of exact solutions of Einstein equations. One 
can even argue that the smooth (de)-phantomization 
is a sufficiently general proper ty of Einstein equation s. 
Therefore if we were (following lAndrianov et aL ( 20051 )) 
guided by a belief that Einstein equations for gravity 
are more fundamental than the concrete form of the 
Lagrangian for other fields, then one can conclude that 
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the "crossing of the phantom division line" is a new 
physically plausible fundamental property of gravita- 
tion. 

The particular interest to models with the phantom 
fields is caused by their prediction of so-called "Cos- 



state, then 

cb{t) = 27r 2 eGDp(t)b 2 (t) 7 



(2) 



mic D oo msday" also known a s big rip [C aldwell et al 
2003h. Istarobinskvl feoool). iFrampton fc Takahashi 



2003 ). lNesseris fc Perivolaropoulosl|2004l) . lNoiiri fc Odint 
2004( ). lGonzalez-Di'aj |2004 ). In case of the phantom 
energy we have w — pj (c 2 p) = — 1 — e with e > 0. Inte- 
gration of the Einstein-Friedmann equation for the flat 
universe results in 



where b(t) is the throat radius of a Morris-Thorne 
wormhole and D is the dimens ionless quantity. Ac- 
cording to iGonzalez-Diaz ( 20041 ) we can choose D ~ 4 
(This is true only if w < 0. If < w < 1 then D ~ A 



a(t) 



p(t) = p Q 



«o 



- (1 _^)2/3e' 

a(t) \ 3e po 
a 



(1) 



(1"€*) S 



where £ = e^/6irGpo. We choose t — as the present 
time, ao ~ 10 28 cm and po to be the present values of 
the scale factor and the density. There, if t = t* = l/£, 
we automatically get the big rip. 

Remark 2. Even though we are certainly uncomfort- 
able with the conclusion that the universe might be 
doomed to be (literally) blown up in the big rip, there 
is simply not too many options available to avert or 
delay it. So far, we only know of four: 

(i) treating phantom energy as nothing more than 
an effective model (see above); 

(ii) using the quantum e ffects to delay the singularity 
( Noiiri fc Odintsovl (|2004l )): 

(iii) introducing a new time variable such that 
the big rip singularity will be delayed indefinitely 
(i.e. its occurrence will be pushed to f -> oo) 
( Gonzalez-Diaz et alj (|2007l )): 

(iv) the last and most intriguing of all: avoiding the 
big rip by means of another comparatively "big" cos- 
molo gical effect - the so-called big trip ( s ee bel ow) . 

In Gonzalez-Diaz fc Jimenez- Madrid! |2004h Pedro 
F. Gonzalez-Diaz, and Jose A. Jime-nez-Madrid had 
introduced a hypothesis that a smooth exit from the 
phantom inflationary phase might actually be achiev- 
able provided one takes into account a multiversity sce- 
nario in which, right before reaching the big rip, the 
primordial phantom universe travels in time towards 
the future state, ultimately landing in a state where it 
will be filled with just usual radiation and no phantom 
fields. This transition had been called "the big trip" 
and it had been assumed it would take place due to a 
sort of anthropic principle, reserving the role of the final 
destination (of the time tra nsition) for our c urren t uni- 
verse. It has been shown in iGonzalez-Diazl (2004) that 



if p = — (1 + e)c 2 p is a cosmological fluid's equation of 



< t Y +3w)( 1+3w '>/ 2w I (4w 3 / 2 ); see lBabichev et al.l d2004l) ). 
The equation (|2|) describes the evolution of b(t) with 
respect to the phantom energy's accretion. Integration 
of the © gives us 



1 2ir 2 ep GDt 



1 

W) ~ b o c(l 

Therefore at 

~ < 
t 



e(cV6irGp + 2vr 2 Po b GD) 



(3) 



(4) 



we get b(t) = oo. As we can see t < t*, and therefore 
this universe will indeed turn achronal well before the 
occurr ence of the big rip. In accord to Gonzalez-Diaz! 
(2004). in the process of phantom energy's accretion, 
for the time period starting right after the radius of the 
wormhole b(t) exceeds the scale factor a(t) and stretch- 
ing up to t, the wormhole becomes an active Einstein- 
Rosen bridge which can, in principle, be used to escape 
the big rip. 

Remark 3. The notion of big t r ip has had its share 
of criticism - the article iFaraonil (|2007l ) of Faraoni is 
by far the most recent and thorough one among those 
expressing the doubts about the model. However, the 
objections of Faraoni has been defended against in a re - 
sponse article Gonzlez-Diaz fc Martm-Morunol ( 2007 ). 
where the detailed answers to all the objections were 
given. It is interesting to note that even phantom mod- 
els with no big rip singularity do admit wormhole so- 
lutions and hence allow for a possib i lity of big trip via 
such a wormhole ( Astashenok et al. ( 20121 )). 

The description of big trip can be executed with help 

of Wheeler-DeWitt equation f or wave function of the 

univer se. To derive it authors of Gonzalez-Diaz fc Jimenez-Mad rid! 
( 20041 ) had introduced the spacetime manifold M for a 
flat FRW universe with metric 



ds 2 = -N 2 dt 2 + a(t) 2 dn 



3< 



(5) 



where N is the lapse function arising due to foliation of 
the manifold M, df2§ is the metric on the unit three- 
sphere, and a(t) is the scale factor. Next, they had 
assumed that the parameter of equation of state w is 
time-dependant, but in such a way that at any moment 
of time we should have 



w = 0, 



(6) 
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so the scalar curvature R would have to be generalized 

by 



R^ R = R- 



aw In a 9 



(7) 



where "dot" indicates differentiation with respect to 
time. Now consider the universe filled with scalar field 
<j) with the Lagrangian 



section we are going to show that by employing the 
isospectral symmetries like the Moutard and Darboux 
transformations, it is possible to construct the exact so- 
lutions of (|13j) . Or, to say it in a more rigorous way, 
there is a simple method of obtaining the "potentials" 
p(u>, a) for which the equation (|13p can be solved ex- 
actly. 



L=\ {d^f - V{<t>) = P = wp. 

The action integral of the manifold M with boundary 
dM will have the form 



S = 



A I 



R 



dxy/=g\-- 



wp 



OAI 



(fxV^hTrK, (8) 



2 Moutard transformation for 
Wheeler-DeWitt equation 

Let us introduce two new independent variables u and 
v. 



u = — + In a, v = — In a. 



(14) 



where K is the extrinsic curvature, g — detg^ and 
h is the determinant of the general threemetric on a 
given hypersurface at the boundary dM. Note, that 
herein we are using the natural system of units where 
8ttG = c = 1. Since yf^g ~ Na 3 one can integrate (JSJ 
over the spatial variables and substitute t — ¥ it, thus 
reducing (|SJ) to Euclidean action: 



and 

w — 2(u + v), a — e~ v . 

As v — > +oo we'll have a — > and v — > — oo would 
imply a — > +oo. The case 



u + v 



1 

2' 



Ndr 



aa 



N 2 



w aa 
2N 2 



wa 3 p 
3 



(9) 



where ' = d/dr. Then in the gauge where N = 1 we 
have Hamiltonian constraint 



51 , s a 3 p 
and the momenta conjugate to a and w: 



1 • 2 
/ | -uia 



lad 



■1. 2 

t-aa 



(10) 



(11) 



However, the rules of quantum description calls for 
them to be replaced by the corresponding quantum op- 
erators: 



7T„ -> I 



d_ 

da ' 



7T„, — > 



.d_ 

dw 



(12) 



By making this change inside of the Hamiltonian con- 
straint we will finally obtain the Wheeler-DeWitt equa- 
tion: 



( d 2 Id 2 3 g i 



(13) 



corresponds to a dS universe and when 

1 

u + v < ——, 

we automatically have a phantom equation of state. 
By direct computation 



d 1 d 



d Id Id 



dw 2 du ' dw 2 da 2 dv 

Therefore, in terms of the new variables the equation 
([T3l has the form 



/ d 2 



\dudv 
where 



+ U(u,v) *(u,w) = 0, 



(15) 



U(u,v) = 3e- 6v p(u,v) 



The equation ([15]) admits the Moutard transformation 
(MT). Namely, let * and $ be two solutions of the (fl5l) : 



1 d 2 ^ 1 <9 2 $ 



= U(u,v). 



^ dudv $ dudv 
Define a 1-form d^f*; $] such that 



(16) 



For all intents and purposes, (1131) is a hard nut 
to crack even for the simplest physically meaningful 
boundary conditions. However, as we shall see, this 
fortress has a hidden back door; in fact, in the next 



d0 $ = — $ - — * - dv( — $ - — * , 
1 cm du J \dv dv J 

(17) 



•1 



and 



= y ef0[tf;$]. 



(18) 



Note that since by definition both \fr and $ are solutions 
of (fT6|) . the one-form is closed, i.e. 

9 2 6>[*;$] _ 2 6»[*;$] 
dudv dvdu 

and thus the shape of the contour of integration V is 
irrelevant. 

The MT has the form 



m _ fl[tt;$] 



?7 (1) = U + 2— — ln$. 
auov 

It means that 

It is worth pointing out at this step that 
$(i) = I 



(19) 



(20) 



(21) 



(22) 



rather then zero. 

The Moutard transformation (|T9j). ([22]) , J20j can be 
iterated several times, and the result can be expressed 
via corresponding Pfafhan forms. 
Example 1. 

Let U = 0, therefore 

* = A(u) + B(v), <f> = a(u) + l3(v), 

with arbitrary functions A(u), B(v), a(u) and (3(v). 
Using (|20|) we get 



2q'(u)fi» 



(23) 



Using ([191) we can calculate the solution of equation 
(|2ip with the potential (1231) . that would contain two 
arbitrary functions and B(v): 

(1) = 2(/3(^)A( M ) - B(v)a(u)) + gi(u) + 9 2 (v) 
a(u)+P(v) 



where 



i(u)= / du(A'(u)a(u) - A(u)a'(u)) , 



(«)= dv(B(v)(3'(v)- B'(v)(3(v)) 



The second solution (without arbitrary functions) is 
given by the formula (f22j) : 



$(!) = 



1 



(24) 



One simple yet interesting example follows if one takes 
a(u) — cosh(w), f3(v) = cosh(w). (25) 



In this case the expression (124)) describes the normaliz- 
able wave function: 



oo r+oo 



dudv\§ (1) \ 2 < 



oo J — oo 



unlike the original solution $, which does not belong 
to L 2 (u,v) for any (non-zero) functions a(u) and j3{v). 
Furthermore, the new solution clearly has a maximum 
at u = v = 0, i.e. at a = 1 and w — (which also 
happens to be a saddle point of potential EA 1 )). Thus, if 
we operate in framework of this model we would have to 
find ourselves with large probability in a dust universe 
with a finite value of scale factor. 

Example 2. Dressing of the dS universe 
Let p = A/3 = const > 0, so U(u,v) = e" 6u A. The 
solution of the (TT5|) has a general form 

*(u,«) = J dnf{K)e KU+ ^ 

with an arbitrary integrable function /(ft). Let us con- 
sider a special case: 



(26) 



with real- valued Kj . For the sake of simplicity, let n = 
2, c(Kj) = 1/2, K2 = —Ki = —k. Then 

A « 

= cosh(7), y(u,v) = ku+ — e , (27) 

OK 

and we can write the result of the Moutard dressing as 
1 



$(i) 



cosh(7) 
and the density p(u, v) as 

,« = A (\ 1 



cosh 7 



(28) 



(29) 



Note, that in this case the maximal probability distri- 
bution will not be localized at one point in the (o, w)- 
phase space, but instead we have a continuous curve 
a(w) that solves the equation j(u, v) = 0. 
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These examples show that, unless we have a very 
special case (like U = 0), the eigenfunctions of the 
Wheeler-DeWitt equation (fTB"]) would in general have 
not just one but many (possibly infinitely many) max- 
imums in the (a, w) space, thus leaving us with the 
burden of choosing the appropriate w that the universe 
should attain after the big trip. Of course, there exist 
different possibilities for such a choice - for example, us- 
ing the final anthropic principle or restricting ourselves 
to a discrete set of possible values of uj (that is treated 
as a quantum- mechanical parameter); nevertheless, it 
is also possible that such a freedom might actually im- 
ply the possible deficiency of the equation itself. In 
fact, as we shall see below, once we take into account 
possibility of w ^ and fc ^ 0, the resulting equation 
would appear to favor one special value of w, namely: 
w = -1/3. 



3 Generalized Wheeler-DeWitt equation and 
the big trip 

The standard definition of a big trip implies that it 
is a cosmological event thought to occur in the future 
during which the entire universe is engulfed inside of a 
gigantic wormhole (hence, the term "big"!), traveling 
through it along space and time. 

However, it would appear that yet another "Big" 
cosmological event might be possible; an event that, by 
analogy, we would like to call a "Big Meeting" . 

In o rder to get to it, we will use the approach 



from iGonzalez-Diaz fc Jimenez-Madridl (|2004l ) with al- 
teration: we would not impose the restrictions w = 
and k = 0. As before, we are assuming that p — p(w, a). 
Differentiating the Einstein-Friedmann equation H 2 = 



(to 



(30) 



where 



R 



6(d 2 + ad + k) 



dp dp 

P = tt« + S~ w - 
da dw 



The corresponding euclidian action is 



I = / drN 



6aa' 3 — aF(a, a! , w, w') 



2a'N 2 



+ 3ka + a 3 wp 



(31) 



where 



F(a, a', w, it/) = 9(1 + w)(a' 3 - ka')- 



Introducing the momenta according to \Y1\ and redefin- 
ing the (classical) Hamiltonian via 7r a , tt w , a and w will 
provide us with a generalized Wheeler-DeWitt equa- 
tion: 

4(a|f+2p+3(l + 3^) 2 ^l- 

(j£) 2 a 6 (l + 3 W )vl/( a ,u;). 



_4a|^^^ = -9 

ow oaow 



(32) 



The WDE (1521 is different from the WDE (U51 . 
obtained in IGonzalez-Diaz fc Jimenez-Madridl (|2004l ) 
since the absence of condition w = has made a fur- 
ther simplification of the actions ©, © not possible. 
As we shall see, it is this new equation (|32|) that results 
in a new effect - the "Big Meeting" . 

First at all, lets consider the special case w = —1/3. 
The right hand side of the (|32|) will be equal to zero. 
Moreover, we will have p = p(a,w = —1/3) ~ a~ 2 , and 
thus 

dp 

a^- + 2p = 0, 
oa 

which reduces the WDE (El to 



logo d 2l t> 
a dadw 



= 0. 



(33) 



Using the power series expansion of function \&(-,a;) 
around the point (•, — 1/3) 



1 f V 
y(a,w) = y(a,-l/3) + ^-c n (a)\w+-) , 

n— 1 ' ^ 



and substituting it into (f33|) would result in dc\ (a) /da 
0, and thus in 



dw 



p/3 — k/a 2 with respect to t will get us a different ex- d^f(a, w) 
pression for the scalar curvature R: 

R = R+^r (9(1 + w)(d 3 + fed) + a 3 p) , 



_i/3 = c\ = const. 



(34) 



On the other hand, equation (1331) is invariant with re- 
spect to the transformation 



*(o,w) *(a,w) - fx(a) - f 2 (w), 



(35) 



for arbitrary (differentiable) functions f\ %. Substitut- 
ing (|33|) into the (f3~4"]l and choosing df2{w)/dw — c\ at 
w = —1/3 we get without loss of generality 



d^(a,w) 
dw 



= -1/3 



0. 



(36) 



Therefore, in the case of general position, for any fixed 
a the function ^(a, w) = <& a (w) has an extremum at 
w = —1/3. Moreover, the condition that function ^> 
has to be normalizable implies that the probability dis- 
tribution I* (a, w)\ 2 for any given value of the scale fac- 
tor has a peak at w = —1/3. In order to prove this 



6 



lets return to equation ([32]) . Consider its solutions that 
satisfy the property 



dw lw=W0 ~ ' 



(37) 



for any given a and some fixed wq = const. Then 
p(a, w ) = p = C 2 a~ 3( - Wo+1 ^ (C = const) and at point 
(•) wo): 



dpo. 



-3(u> + l)po, 



a — u)=«)o = -31oga/9o- 



Besides 
<9 2 * , 



dadw lw=W0 ~ da { dw lw=W0 



(38) 



(39) 



(40) 



Substituting 



([3"9"f and gDJl into the (02} we get 



(f + 3w Q )(Aa~^+^(l - 3fcC- 2 a 1+3 ^)0U, =u;o 
-81 C 2 a- 6w ° log 2 a*) = 0. 

(41) 

Using (|4Tj) one can conclude that either w = — 1/3 or 
we have to satisfy the following equation 



(l-3fcC- 2 a 1 + 3 «'»)0| 



81C" 



a -3(w -l) 1q 2 a4f _ 



(42) 



If = the general solution of the (|4"2"]l has the form 

V(w ,a) = ci/o(^) + c 2 K (z), (43) 

where Ci 2 are arbitrary constants, z = SCa 3 ' 1- " 10 ^ 2 , 7 
and Ko are the Bessel functions of the first and second 
kind correspondingly. This function has to be normal- 
izable: 



da\^(w ,a)\ 2 < +00, 



(44) 



which implies that c\ — 0. Substituting (143)) into the 
(|40|) results in 

a* I _ 

dw \ ™=1>>0 ~ 

= -^a^-^/HogaKt (ZCa^ 1 -^' 2 ) = 0, 

which can only be satisfied for arbitrary a if and only 
if wq — 1 and C is the solution of the equation 

CKi(3C) = 0. 

But if wo = 1 then the wave function '5 = const (and 
the same will be the case for the density p (see (f3"2"j) )). 



But the only normalizablc function of this type is ^ = 
0. Therefore, the only possible non-trivial solution that 
comply with (|37j) in framework of normalizable wave 
functions of the universe correspond to the case wo — 
-1/3. 

The end of the proof. 

Now lets say a few words on account of possible ex- 
istence of other peaks M* = (a*, it;*) on the graph of 
solutions of the (|32|) that satisfy 



5* 9* 



and 



A 



9 2 *a 2 * 



dw 2 da 2 \ dwda J 



> 0. 



(45) 



Surprisingly, it is possible to eradicate these extrema. 
To show it, one should take into account the dimen- 
sionless variables in expressions like log a. Keeping in 
mind that our universe was born via big trip from the 
paternal universe with scale factor L, we have to re- 
place log a — > log(a/L). On the other hand, we expect 
that right after the big trip the initial value of the scale 
factor dj will be a, ~ L. Lets consider such a universe 
with w = const. In order to estimate the probability 
to find this universe just after the big trip (that is with 
a t ~ L ~ a*) one must use the WDE (J32J which reduces 
to a simple form 

-4a- 8 <™ +1 >(3«, + 1)01^=0. 
If w j= -1/3 then (see (05])) 



A = 



/ dH y 

V dwda I 



\M, 



<o, 



and we have no extremum at all. 

Note one interesting feature of a universe with w — 
— 1/3. If w = — 1/3 the solution for scale factor has a 
simple form 

a(t) = At + a , 

where A and ao are constants. For the expanding uni- 
verse A > 0. It is possible to join the metric with this 
scale factor to any FRW universe at any time in the 
expanding phase. Consider for example a flat de Sitter 
universe: 

a(t) = Cexp(Ai) 



where C — const and A = y^SirGpA/S. The matching 
conditions for arbitrary t = r require continuity of the 
scale factor together with its first derivatives: 
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ci(t) = C exp(Ar) = At + ao, 

d(r) = CAexp(Ar) = A. 

For ao we have clq = C(l — Ar) exp(Ar). Therefore the 
joining is trivial: the continuity condition of continu- 
ity of a(t) uniquely determines the parameter A and 
the continuity condition of a(t) does the same for the 
constant ao- 

4 Conclusion 

The simple method based on Moutard transformation 
for constructing exact solutions of the Wheeler-DcWitt 
equation (fl~3)) is presented. A direct application of this 
method has demonstrated that in general, unless some 
kind of additional restriction on the possible values of w 
is imposed, we would have to face a problem of somehow 
choosing a suitable w from (in general) infinite set of 
possibilities whenever we try to describe the universe 
after the big trip. 

However, it appears that should one consider the 
generalized form of Wheeler-DeWitt equation (|32|) in- 
stead, the conclusion would be drastically different. In 
this case the probability distribution has the peak(s) 
at w = —1/3 for any a provided that we consider the 
universe right after big trip. By analogy with the "big 
trip" this result can be called a "Big Meeting" . What 
it means is that in multiverse the vast majority of uni- 
verses experiencing the big trip should end up in a state 
with w = —1/3 right afterwards. 
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